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pumping 

A. W. Carr and M. Saffman^ 

^Department of Physics, 1150 University Avenue, 

University of Wisconsin, Madison, Wisconsin 53706 

(Dated: April 24, 2013) 

We propose and analyze an approach for preparation of high fidelity entanglement and antiferro- 
magnetic states using Rydberg mediated interactions with dissipation. Using asymmetric Rydberg 
interactions the two-atom Bell singlet is a dark state of the Rydberg pumping process. Master 
equation simulations demonstrate Bell singlet preparation fidelity T = 0.998. Preparation of anti- 
ferromagnetic spin plaquettes is also demonstrated. 
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Neutral atoms are providing a new tool for studying 
many body physics and quantum magnetism [1|. Re- 
cent experiments have probed strongly correlated spin 
systems relying on short range contact interactions be- 
tween cold atoms |2|] , or Coulomb interactions of trapped 
ions[3|. Much recent interest has focused on antiferro- 
magnetically ordered spin states. When the spin inter- 
actions are mediated by short-range scattering it is chal- 
lenging to reach the extremely low temperatures needed 
to observe antiferromagnetic ordering y. Recent progress 
with Rydberg excited atoms [sj has demonstrated strong 
and long-range dipolar interactions which are suitable for 
creating magnetic phases with long range order [6|. 

In this letter we propose and analyze an approach to 
entanglement generation and spin ordering which relies 
on dissipative dynamics with Rydberg state mediated in- 
teractions. It is well known that dissipative dynamics can 
be used for creating entanglement [7], and more gener- 
ally for universal quantum computational tas ks fei] as has 
been demonstrated in recent experiments |9[ llQj. These 
developments have led to a high level of activity on this 
topic resulting in approaches to dissipative entanglement 
generation in a range of physical settings [11 1. 



A dissipative approach to antiferromagnetic ordering 
using Rydberg interactions was proposed in[l2|. Here we 
consider an arrangement with spin-dependent Rydberg 
interactions, which allows us to prepare the two- atom 
spin singlet state as a dark state of the dissipative evo- 
lution. In contrast to coherent blockade experiments [13| 
which rely on minimizing dissipation in order to maxi- 
mize the fidelity of the target quantum state, the present 
approach exploits spontaneous emission, yet can be used 
to prepare a maximally entangled singlet state with fi- 
delity exceeding 0.998. This dynamics enables high fi- 
delity entanglement at long range which will be useful 
for teleporting gates in a spatially extended qubit array, 
as well as creating strongly correlated spin systems. Re- 
markably the dissipative approach described here is ca- 
pable of creating the same entanglement fidelity as the 
coherent Rydberg blockade gate[l^ but with 1500 times 
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FIG. 1. Level diagram for pumping into an asymmetric 
ground state using unequal Rydberg interactions. The boxed 
inset shows how tt polarized light can be used to recycle all 
other Zeeman ground states in the example of Cs atoms. See 
text for details. 



smaller Rydberg interaction. This implies that entangle- 
ment can be extended to much longer interparticle sepa- 
rations which will enable efficient computation and many 
particle entanglement in extended qubit arrays. 

Our approach to dissipative preparation of a spin sin- 
glet is illustrated in Fig. [TJ Consider two atoms with 
ground states |1), |2) and Rydberg states |3), |4). State 1 
is coupled to 3 with Rabi frequency fti and 2 is coupled 
to 4 with 1^2- The driving fields are detuned by A. The 
energy splitting of the ground states hw2i is assumed to 
be a few GHz corresponding to atomic hyperfine ground 
states. The splitting is sufficiently large for us to neglect 
off-resonant coupling of |1) -^ |4) and |2) -^ |3). States 
3, 4 decay to 1, 2 with rate 7 and equal branching ratios. 
In addition there is a field coupling |1) ^ |2) with Rabi 
frequency ftg. This could be implemented optically with 



two-photon Raman transitions or with a microwave field. the Schrodinger equation takes the form 



The Rydberg states interact with energy V33 = V44 
when both atoms are excited to the same state and en- 
ergy V34 = V43 when the atoms are excited to differ- 
ent states. In general the Rydberg interaction energy 
depends on angular degrees of freedom so the coupling 
strength will be different when |3), |4) correspond to dif- 
ferent Zeeman sublevels|15|. The idea of using asym- 
metric Rydberg couplings for quantum state control was 
previously used in several papers [16-18], In those pa- 
pers a very strong interaction asymmetry was required 
for high fidelity control. The idea described here re- 
quires only a weak asymmetry to enable a multi-atom 
optical pumping process. We now choose a detuning 
A = ^3/2fi = A33/2. Thus when the ground atomic 
state is |33) or |44) both atoms will be excited and spon- 
taneously decay from the Rydberg level, thereby popu- 
lating ah four ground states |11), |12), |21), |22). On the 
other hand excitation of |12) or |21) will be off-resonant 
by an amount S = A33 — A34, and excitation of a single 
atom will be off-resonant by A. Provided the off-resonant 
excitation rates are small compared to excitation of the 
symmetric ground states the system will be pumped into 
an asymmetric two-atom state. 



Although spontaneous emission in real atoms popu- 
lates other hyperfine ground states not in the basis |1), |2) 
this can be dealt with using recycling lasers as shown in 
the inset of Fig. [H Explicitly for the example of Cs atoms 
we take |1) = |/,m) = |3,0), |2) = |4,0) and add tt po- 
larized lasers coupling |65i/2,/ = 3) ^ |6pi/2,/ = 3), 
I6S1/27/ = 4) ^ |6pi/2,/ = 4). These lasers do not 
disturb |1), |2) but recycle all other ground states. 



We can gain some insight into appropriate parameters 
by considering the limit of coherent Schrodinger evolu- 
tion. Assume two, two-level atoms with ground state 
|1), Rydberg state |3) and Rydberg interaction strength 
V33 = hAss. There is resonant excitation of the dou- 
bly occupied Rydberg state when the laser detuning is 
A = A33/2. The states |13) and |31) are off-resonant by 
A33/2 and are only weakly excited. Writing the state 
vector as 
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Here s = ^(ci3 
and A = uj 



C31), uja is the transition frequency, 
uja is the laser detuning from the non- 
interacting atomic transition resonance at uja- Adia- 
batically eliminating the singly excited amplitude s at 
A = A33/2 we get an effective two- level system for states 
1 11) and 1 33) resonantly coupled with the Rabi frequency 
Qr = {\^Q)'^/Ass. Allowing for four- levels as in Fig. [1] 
the symmetric states are resonantly coupled to Rydberg 
levels with Qr while the antisymmetric states are also 
coupled with Qr^ but at a detuning S = A33 — A34. 

In the absence of symmetry breaking we expect the 
atoms to end up in a state 1?/^) ~ |12) + e*'^|21) with 
(p an undetermined phase. In order to create coherence 
between |12) and |21) we add the transverse drive Qg. 
The singlet state |— ) = 4=(|12) — |21)) is invariant with 
respect to X rotations, whereas the m = triplet state 
1+) = "7^(112) + |21)) couples to states with m = ±1 
which are subject to repumping via the Rydberg states. 
We thus expect the combined action of Qi^Q2^^g with 
the detunings and Rydberg couplings specified above will 
drive the atoms into the maximally entangled "antiferro- 
magnetic" state |— ). 

Denoting the probability for the atoms to be in the 
antiferromagnetic state |— ) by Paf the rates at which 
probability enters and leaves this state due to one and 
two-atom excitation processes are 



ri,in = (1 - Paf)2 



n,out = Paf2- 
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We then solve n^ia + ^2,in = ''i.out + ?'2,out, to find the 
equilibrium population 
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FIG. 2. (color online) Comparison of Paf from Eq. ([2]) 
(solid line) with singlet state fidelity FbgII from numerical so- 
lution of Eq. Q (red dots). Parameters: Q/27r = 0.01 MHz, 
7 = l./(0.73 ms), A34/A33 = 0.2, Qr = 2^2/^33, and 
Qg — O.Sr^i?. The atomic parameters correspond to the Cs 
125pi/2 state with |3), |4) rrij — ±1/2 and the quantization 
axis perpendicular to the line containing the atoms. The 
inset shows the time dependence for A33/27r = 1 MHz. 



This expression is plotted in Fig. [2] for representative 
parameters. As the interaction strength increases the 
fidelity of the Bell singlet tends to unity. This is also 
borne out by the limiting expression found by taking 
^ ::^ 7, Vtn^ A33 ^ Vt so that both the one and two 
atom rates out of the singlet state are small, which leads 

to Paf ^ (1 + 2f^) / (1 + IT^) • We see that provided 
1^ > 7 then Paf ^ 1- 



To verify the approximate rate equation analysis we 
have numerically solved the two-atom master equation 



t->-\-'- 



(3) 



with H = Hi (g) /2 + /i H2 + V, £ = £1 /2 + /i /:2, 

and /i, I2 are 4x4 identity matrices. Working in the in- 
teraction picture and making the rotating wave approx- 
imation the one atom operators expressed in the basis 

{|l),|2>,|3>,|4>}are 
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Numerical solutions of ([3j) are used to extract the fidelity 
of the Bell singlet state | — ) given by F = ^(pii;22 + 
P22;ii) + |pi2;2i|- Valucs of F fouud from integrating 
to t = 90 /Qr starting from the initial condition |11) 
are compared with the approximate result for Paf that 
comes from solving Eqs. ([T]) in Fig. [21 

We see that the approximate result agrees well with 
numerical solutions. The maximum Bell state fidelity is 
i^Beii = 0.9988 at A33/27r = 3 MHz. This is essentially 
the same fidelity as the best found in [14] for the coherent 
Rydberg blockade controlled phase gate. It is notewor- 
thy that the dissipative approach does not require single 
atom addressing and the fidelity is achieved with about 
1500 times smaller Rydberg interaction strength. This 
implies that high fidelity entanglement can be achieved 
at very much larger atomic separations than for the co- 
herent interaction. A viable approach to long range 
gates in an array of qubits could thus be based on 
teleportation[19| using the dissipative mechanism for es- 
tablishing entanglement, followed by short range coher- 
ent gates between neighboring qubits. 

Our calculations ignore undesired entanglement be- 
tween spin and center of mass degrees of freedom. This 
can be suppressed, despite the presence of spontaneous 
emission from the Rydberg levels, provided we confine 
the atoms in the Lamb-Dicke regime and use magic 
ground- Rydberg trapping potentials [20]. 

The antiferromagnetic state with neighboring spins an- 
tiparallel is trivially the ground state of the Ising model 
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operators acting on atom j, 731 = 732 = 741 = 742 = 7/2 
and all other ^yki = 0. The interaction term is 

V = hAs3 (|33)(33| + |44)(44|) + hAs4 (|34)(34| + |43)(43|) 



when we restrict to nearest neighbor couplings and there 
is no frustration. Finding the ground state of H on a 2D 
lattice with couplings that extend beyond nearest neigh- 
bors and with a local transverse field (corresponding to 
Qg) is generally a hard computational problem [2l| which 
may be amenable to simulation using the Rydberg cou- 
plings described above. 

Numerical solution of the master equation for more 
than several atoms requires large computational re- 
sources. We have instead used a Monte Carlo wavefunc- 
tion technique [2 2| to study the dynamics on three and 
four spin plaquettes. We will refer to internal states |1), 
1 2) as spin up and down. In the case of three spins on 
an equilateral triangle the dynamics is frustrated and we 
find equal superpositions of the six distinct states with 
M = ±1/2. Results for four atoms on a square plaquette 
with rtg = are given in Fig. [3l In this case we pump 
close to 100% of the population into a superposition of 
the M = states \a) = |1212) and \b) = |2121), corre- 
sponding to the ground state of the Ising model (J5]). The 
remaining few percent of the population is fairly evenly 
distributed over the 12 higher energy states. Solution 
of ([5]) on a square plaquette with 1/R^ van der Waals 
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FIG. 3. Population of AF± states from Monte- Carlo sim- 
ulations of Eqs. (|3l4p on a square plaquette and averaged 
over 100 trajectories. The blue and red curves give the pop- 
ulations of 1 1212) and |2121). Numerical parameters were 
Q/27r = 0.02 MHz, 7 = l./(0.3 ms), A33/27r = 0.4 MHz, 
A34/A33 = 0.85, A = A33/2, Qr = 2QVA33, and Qg = 0. 
The coupling strengths between opposite corners were A33 = 
A33/8 and A^4 = A34/8. 



scaling of the interaction strength shows that the ground 
states have energy ~y^ , while the next four states are at 
^. Our physical model with dissipation and Rydberg 
interactions shows relaxation to the ground state with 
purity of about 96%. 

The calculation was done with Qg = so no coherence 
is generated between |a), |6). Setting Qg ^ gives non- 
stationary dynamics since the J = M = state of four 
spins cannot be expressed as a superposition of states 
I a), \b) for which neighboring spins are antiparallel, and 
global rotations do not lead to a stationary dark state. 

In summary we have shown how the combination of 
weak Rydberg interactions and global rotations pump 
pairs of atoms into highly entangled singlet states. A re- 
lated approach to entanglement generation using electro- 
magnetically induced transparency with dissipative dy- 
namics has been presented in [23|. The entangled sin- 
glet states can then be used as a resource for high fi- 
delity, long distance gates in a qubit array based on gate 
teleportation|19|. Going beyond pairs of atoms we show 
that the dynamics mainly pumps atoms into an inco- 
herent superposition of spin zero states corresponding to 
the ground state of an antiferromagnetic Ising interac- 
tion. The overlap of the multi-atom states so generated 
with the ground states of long range Ising models and 
their dynamical evolution are open questions for future 
studies. 
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